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We consider the combined influence of linear damping and noise on a dynamical finite-time- 
singularity model for a single degree of freedom. We find that the noise effectively resolves the 
finite-time-singularity and replaces it by a first-passage-time or absorbing state distribution with a 
peak at the singularity and a long time tail. The damping introduces a characteristic cross-over 
time. In the early time regime the probability distribution and first-passage-time distribution show a 
power law behavior with scaling exponent depending on the ratio of the non linear coupling strength 
to the noise strength. In the late time regime the behavior is controlled by the damping. The study 
might be of relevance in the context of hydrodynamics on a nanometer scale, in material physics, 
and in biophysics. 
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I. INTRODUCTION 

The influence of noise on the behavior of nonlinear dynamical system is a recurrent theme in modern statistical 
physics 0. In a particular class of systems the nonlinear character gives rise to finite-time-singularities, that is 
solutions which cease to be valid beyond a particular finite time span. One encounters finite-time-singularities in 
stellar structure, turbulent flow, and bacterial growth 

H H H The phenomenon is also seen in Euler flows and in 
free-surface-flows 0, @. Finally, finite-time-singularities are encountered in modelling in econophysics, geophysics, 
and material physics |7l H H Ha ITU . 

In the context of hydrodynamical flow on a nanoscale |l2j| , where microscopic degrees of freedom come into play, it 
is a relevant issue how noise influences the hydrodynamical behavior near a finite-time-singularity. Leaving aside the 
issue of the detailed reduction of the hydrodynamical equations to a nanoscale and the influence of noise on this scale 
to further study, we assume in the present context that a single variable or "reaction coordinate" effectively captures 
the interplay between the singularity and the noise. 

In a recent paper jl3l | we investigated a simple generic model system with one degree of freedom governed by a 
nonlinear Langevin equation driven by Gaussian white noise, 

dr \ 

*=-2j^ + 1> <*>>(') = **(*)• (I'D 

The model is characterized by the coupling parameter A, determining the amplitude of the singular term, the index 
fi > 0, characterizing the nature of the singularity, and the noise parameter A, determining the strength of the noise 
correlations. Specifically, in the case of a thermal environment at temperature T the noise strength A oc T. 

In the absence of noise this model exhibits a finite-time-singularity at a time to, where the variable x vanishes 
with a power law behavior determined by /i. When noise is added the finite-time-singulari ty e vent at to becomes a 
statistical event and is conveniently characterized by a first-passage-time distribution W(t) For vanishing noise 
we have Wit) = Sit — to), restating the presence of the finite-time-singularity. In the presence of noise Wit) develops 
a peak about t = to, vanishes at short times, and acquires a long time tail. 

The model in Eq. 1)1.1(1 has also been studied in the context of persistence distributions related to the nonequilibrium 
critical dynamics of the two-dimensional XY model 01 an d in the context of non-Gaussian Markov processes (lt| . 
Finally, regularized for small x, the model enters in connection with an analysis of long-range correlated stationary 
processes |l7| . 
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From our analysis in ref. [l3j it followed that for n = 0, the logarithmic case, the distribution at long times is given 
by the power law behavior 

W(t)~t~ a , «=^ + ^- (1.2) 

For vanishing nonlinearity, i.e., A = 0, the finite-time-singularity is absent and the Langevin equation l|l.l|l describes 
a simple random walk of the reaction coordinate, yielding the well-known exponent a = 3/2 [iH Il8l llflj. In the 
nonlinear case with a finite-time-singularity the exponent attains a non-universal correction, depending on the ratio 
of the nonlinear strength to the strength of the noise; for a thermal environment the correction is proportional to 1/T. 
In the generic case for /i > we found that the fall-off is slower and that the correction to the random walk result is 
given by a stretched exponential 

W(t) ~ t- 3/2 exp[~A(t-» /(2+ ^ - 1)], (1.3) 

where A — ► A/A/i for [i— * 0. 

In our studies so far we have ignored damping. It is, however, clear that in realistic physical situations friction or 
damping must enter on the same footing as the noise. This follows from the Einstein relation or more generally from 
the fluctuation-dissipation theorem relating the damping to the noise. In the present paper we attempt to amend 
this situation and thus proceed to extend the analysis in ref. 01 to the case of linear damping. We shall here only 
consider the logarithmic case for fi = 0. 

For this purpose we consider the following model for one degree of freedom: 

^ = - T e--^ + ?7, ( W )=A5(i). (1.4) 
at 2\x\ 

In addition to the coupling parameter A, and the noise parameter A, this model is also characterized by the damping 
constant 7. Assuming for convenience a dimcnsionless variable x the coupling and the noise strengths A and A have 
the dimension 1/time. The ratios A/ A and 7/ A are thus dimensionless parameters characterizing the behavior of the 
system. 

It follows from our analysis below that the damping constant sets an inverse time scale I/7. At intermediate time 
scales for jt <C 1 the distribution exhibits the same power law behavior as in the undamped case given by Eq. (|1.2fl . 
At long times for ^t ^> 1 the distribution falls off exponentially with time constant 1/7(1 + A/A), i.e., 

W(t) oc exp[-7(l + A/A)t]. (1.5) 

The paper is organized in the following manner. In Section [D] we introduce the finite-time-singularity model with 
linear damping and discuss its properties. In Section IIHI we review the weak noise WKB phase space approach 
to the Fokker-Planck equation, apply it to the finite-time-singularity problem with damping, discuss the associated 
dynamical phase space problem and the long time properties of the distributions. In Section llVI we derive an exact 
solution of the Fokker-Planck equation and present an expression for the first-passage-time distribution. In Section 
[V]we present a summary and a conclusion. In the present treatment we draw heavily on the analysis in ref. 0|. In 
Appendix [S] aspects of the exact solution are discussed in more detail; in Appendix [5] we consider the weak noise 
limit of the exact solution. 



II. MODEL 



In terms of a free energy or potential F we can express Eq. (|1.4f) in the form 

dx 1 dF ^ . ^ 
dt 2 dx ' 

where F has the form 

F = 7a; 2 + Aln|4 (2.2) 

The free energy has a logarithmic sink and drives x to the absorbing state x = 0. For large x the free energy has 
the form of a harmonic well potential confining the motion. In Fig. have depicted the free energy in the various 
cases. 
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A. The noiseless case 



In the case of vanishing noise Eq. (|2.1() is readily solved. We obtain 



x{t) 



2 7 



-,1/2 



with a finite-time-singularity at 



to 



2l 



log 



,2 7 (t -t) 



1+ x x ° 



1/2 



(2.3) 



(2.4) 



The initial value of x is xq at time t = 0. In the presence of damping a; initially falls off exponentially due to the 
confining harmonic potential with a time constant I/7. For times beyond I/7 the nonlinear term takes over and 
drives x to zero at time to, i.e., x falls into the sink in F. In Fig.^we have shown the noiseless solution x(t). 



B. The noisy case 

Summarizing the discussion in ref. |13| . the stochastic aspects of the finite-time-singularity in the presence of 
noise are analyzed by focusing on the time dependent probability distribution P(x, t) and the derived first-passage- 
time or absorbing state probability distribution W(t). The distribution P(x,t) is defined according to |19t |20| 
P(y,t) = (5(y — x(t))) where x is a stochastic solution of Eq. 1|2.I|) and (■ • •) indicates an average over the noise r\ 
driving x. In the absence of noise P(y,t) — 8{y — x(t)), where x is the deterministic solution given by Eq. I|2.3I) and 
depicted in Fig^ At time t = the variable x evolves from the initial condition xo, implying the boundary condition 
P(x,0) =5{x-xq). 

At short times x is close to xo and the singular term and the damping term are not yet operational. In this regime 
we obtain ordinary random walk with the Gaussian distribution, P{x,t) = (2TrAt)~ 1 ^ 2 exp [— (x — xq) 2 /2At] . At a 
time scale given by I/7 the damping drives x towards a stationary distribution given by P oc exp(— F/A). However, 
at longer times beyond the scale I/7, the barrier X/2x comes into play preventing x from crossing the absorbing state 
x = 0. This is, however, a random event which can occur at an arbitrary time instant, i.e., the finite-time-singularity 
at tg in the deterministic case is effectively resolved in the noisy case. For not too large noise strength the distribution 
is peaked about the noiseless solution and vanishes for x — ► 0, corresponding to the absorbing state, implying the 
boundary condition 



P(0,t) = 



(2.5) 



In order to model a possible experimental situation the first-passage-time or absorbing state distribution W(t) is of 
more direct interest [2(J, |2l| ■ 

Since P(0, t) — for all t due to the absorbing state, the probability that x is not reaching x — in time t is thus given 
by J P(x, i)dx, implying that the probability —dW that x does reach x = in time t is —dW = — J dxdt(dP/ dt) , 
yielding the absorbing state distribution W(t) — — J Q °° dx dP(x,t)/dt [l^. In the absence of noise P(x, t) = S(x—x(t)) 
and W(t) = S(t — to), in accordance with the finite time singularity at t = to. For weak noise W(t) peaks about to 
with vanishing tails for small t and large t. 

The distribution P(x, t) satisfies the Fokker-Planck equation [20ll2l| 



dP 



ld_ 

2dx 



—P + A— 

dx dx 



(2.6) 



in the present case subject to the boundary conditions P(x,0) = 6(x — xq) and P(0,t) = 0. The Fokker-Planck 
equation has the form of a conservation law dP/dt+ d J/dx = 0, defining the probability current J = (l/2)(dF/dx)P— 
(l/2)AdP/dx and we obtain for W(t) the expression 



W(t) = 



dF ' A <9P 

dx dx 



(2.7) 



to be used in our further analysis. Note that there is a sign error in Eq. (3.8) in ref. |13|. 
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III. WEAK NOISE APPROACH 

In this section we apply a weak noise canonical phase space approach to the damped finite-time-singularity model 
and infer the general long time behavior. The treatment follows closely the analysis in ref. |13| . 

A. The phase space method 

From a structural point of view the Fokker-Planck equation H2.6J) has the form of an imaginary-time Schrodinger 
equation AdP/dt = HP, driven by the Hamiltonian or Liouvillian H. The noise strength A plays the role of an 
effective Planck constant and P corresponds to the wavefunction. Drawing on this parallel we have in recent work in 
the context of the Kardar-Parisi-Zhang equation for a growing interface elaborated on a weak noise nonperturbative 
WKB phase space approach to a generic Fokker-Planck equation for extended system [2^, [53, 13 • I n the case of 
a single degree of freedom this method amounts to the eikonal approximation |19l l2lL l25l , see also |26l 1271 ] . For 
systems with many degrees of freedom the method has for example been exp ounded in [28| , based on the functional 
formulation of the Langevin equation |2^,|33|. In the present formulation |22t I23I |24| the emphasis is on the canonical 
phase space analysis and the use of dynamical system theory [^J • 

The weak noise WKB approximation corresponds to the ansatz P cx exp[— S/ A}. The weight function or action S 
then to leading asymptotic order in A satisfies a Hamilton- Jacobi equation dS/dt + H = which in turn implies a 
principle of least action and Hamiltonian equations of motion |33, y4| . In the present context the Hamiltonian takes 
the form 

H=\p{p---2 1X l (3.1) 
2 x 

yielding the Hamilton equations of motion 

(3.2) 
(3.3) 

These equations replace the Langevin equation (|1.4f> with the noise r\ represented by the momentum p = dS/ dx, 
conjugate to x. The equations (|3.2|) and l|3.3|l determine orbits in a canonical phase space spanned by x and p. Since 
the system is conserved the orbits lie on the constant energy manifold(s) given by E = H. The action associated with 
an orbit from xq to x in time t has the form 



dx 




A 


lit 


= —721 — 


Tx +V ' 


dp 


1 


A 


~dt 


= 7P~2 


x 2 ^ 



S(xq —> x,t) = I dt 
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dx 



According to the ansatz the probability distribution is then given by 

S(xq —>■ x, i) 



P(x, t) — P(xq — ► x, t) cx exp 



(3.4) 



(3.5) 



B. Long time orbits 

The zero-energy manifolds delimiting the phase space orbits follow from Eq. I|3.1|) and are given p — and 
p = 2jx + X/2x. The p = sub-manifold corresponds to the noiseless or deterministic case discussed above. The 
p = X/x + 27a; sub-manifold corresponds to the noisy case. By insertion in Eq. I|3.3|) we obtain dx/dt — jx + A/2x, 
i.e., the motion on the noisy sub-manifold is time reversed of the motion on the noiseless sub- manifold. The orbit 
structure in phase space is moreover controlled by the hyperbolic fixed point at (x*,p*) = [(A/27) 1 / 2 , (27A) 1 / 2 ]. The 
heteroclinic orbits passing through the fixed point are given by p — X/x and p — 2~/x and the energy of the invariant 
manifold is E* = — 7A. In Fig [2 we have depicted the phase space with the zero-energy manifolds, the fixed point, 
the heteroclinic orbits and some characteristic orbits. 

The long time behavior of the distribution is determined by an orbit from xq to x traversed in time t. In the long 
time limit this orbit must pass close to the hyperbolic fixed point. Note that in the limit 7^0 the fixed point 
migrates to infinity in the x direction and the long time orbits approach the zero energy sub-manifolds which thus 
determine the asymptotic properties as discussed in ref. |l3[ . 
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Independent of whether the initial value xq is greater or smaller than the fixed point value x* , the long time orbit 
follows the invariant p = X/x manifold towards the fixed point. At the fixed point the orbit slows down and then 
speeds up again as the orbit follows the other invariant manifold p = 2jx towards the endpoint x reached in time t. 
This behavior is also depicted in Fig El 

This scenario allows a simple analysis of the long time behavior of the distribution P(x, t) and the first-passage- 
distribution W(t). Close to the invariant manifolds with energy E* = — 7 A the action associated with an orbit from 
xq to x follows from Eq. (|3.4|) and is given by 



S = -EH+ [ dx- + I d.r-2'.r. 

Jxo x 

or, denoting the relevant manifolds by a subscripts, see Fig. [21 



S = \-ft + A log 



1 / 2 *2\ 
+ J(X -X jll. 



(3.6) 



(3.7) 



At long times we only have to consider the contribution from the orbit leading up to the fixed point. Inserting the 
manifold condition p = X/x in the equation of motion (|3.2[) we thus obtain dx/ dt = —72; + X/2x with solution 



x{tf 



cle^ + x^il-e' 2 ^). 



(3.8) 



C. Discussion 



It follows from Eq. (|3.8|l that the damping 7 sets an inverse time scale delimiting two kinds of characteristic 
behavior. First, for t — > 00 the orbit approaches the fixed point x*. For 7< > 1 we have x 2 = x* 2 (l — exp(— 2jt)) 
and x approaches the fixed point in an exponential fashion. On the other hand, in the intermediate time region for 
jt <C 1 and for jt ^> Xq and Xt 3> x 2 we obtain x 2 = Xq + 2tjx* = x\ + Xt ~ Xt. 

By insertion in the expression l|3.7|l for the action we then obtain in the late time regime for jt S> 1 

5(t)~A 7 i-^, (3.9) 

yielding the distribution and ensuing first-passage-time distribution 

P(t) cx W(t) cx exp(-A7*/A). (3.10) 
Likewise, we have in the intermediate time regime jt <C 1 

S(i)~A7*+~log|i|, (3.11) 
giving rise to the distribution and first-passage-time distribution 

P(t) cx W(t) cx \t\-&. (3.12) 

These results hold in the weak noise limit. We note that at long times W(t) falls of exponentially with a time 
constant given by A/A7. In the intermediate time regime W{t) exhibits a power law behavior with exponent — A/2A, 
independent of I/7 defining the cross-over time. These results will also be recovered from the exact solution discussed 
in the next sections. 



IV. EXACT SOLUTION 



In this section we return to the Fokker-Planck equation l|2.6|l and present an exact solution. This solution is an 
extension of the solution presented in ref. [Iflj and the analysis proceeds in much the same way. Details are discussed 
in Appendices lAl and iBl 
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A. Quantum particle in a harmonic potential with centrifugal barrier 



The Fokker Planck equation has the form 

dP A d 2 P 



7X ■ 



A 

2x 



dP 

dx 



+ 7- 



2 dx 2 

Eliminating the first order term by means of the gauge transformation 

exp(h) = |x|- A / 2A e-^ 2 / 2A , 

we can express the equation in the form 

d 



A 

2X 2 



P. 



A- [exp(- 



-h)P] = H [exp(-/i)P] 



where the Hamiltonian H driving P is given by 

A 2 



H 



1 A2 9 2 
-A 

2 dx 2 



1 



2A 

X 



A7 



1 



7 2 2 

X . 

2 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



This Hamiltonian describes the motion of a unit mass quantum particle in one dimension in a harmonic potential 
subject to a centrifugal barrier of stren gth (A 2 /8)(l + 2A/A) at the origin; A plays the role of an effective Planck 
constant. Note that in Eq. (6.4) in ref. the factor A/2 should read 1/A. 

For A = and 7 = both the barrier and the confining potential are absent; the spectrum of H forms a band and 
the particle can move over the whole axis. This case corresponds to ordinary random walk [19|. Incorporating the 
absorbing state condition in Eq. (|2.5I) by means of the method of mirrors we obtain the results presented in ref. |13| , 
i.e., 



P(x,t) = (2nAt) 



-1/2 





(x-x ) 2 ' 




(x + x ) 2 ' 




exp 


2At 


— exp 


2At 





in the half space x > 0, and for the absorbing state distribution 

• 2 x 1/2 

7T 



W(t) 



x Q exp(-xl/2At)(At)' 3/2 . 



(4.5) 



(4.6) 



For A ^ and 7 = the particle cannot cross the barrier and is confined to either half space; this corresponds to 
the case of a finite-time-singularity subject to noise and an absorbing state at x = and was discussed in detail in 
ref. 0; for reference we give the obtained results below. Note that x and xo should be interchanged in Eq. (6.5) 
and that a factor A is missing in Eq. (6.6) in ref. fl3| . 



P(x,t) 
W(t) 



A/2A+1/2 exp (- 

V 



2Ai 



2-A/2A-1/2 

2A4 +A/A 



j /xx \ 
H2S V A< / 



At 

■exp(-x 2 /2Ai)(2A£) 



(4.7) 



(4.8) 



r(l/2 + A/2A) 

In the present case for A 7^ and 7 7^ the problem corresponds to the motion of a particle in a harmonic potential 
with a centrifugal barrier at x = 0. The spectrum is discrete and becomes continuous for 7 = 0. The pro blem 
is readily analyzed in terms of confluent hypergeometric functions, more specifically Laguerre polynomials |35l l36j| . 
Incorporating the initial condition P(x,0) = 5{x — xq) and introducing the time scaled variables 



x = x, exp(— r yt/2), 
xq = x cxp(+7i/ 2 ), 



we find for P(x,t) 



P(x,t) 



-A/2A+1/2 



7e 



■yt/2 



■ exp 



JA/2A-1/2 A s inh 7 i 
and correspondingly for the absorbing state distribution 



l(x 2 + go) 
2Asinh7t 



7 xx 



2 + 2A y A sinh7< 



W(t) 



2Ax 



-1+A/A 







r(l/2 + A/2A) 



exp 



7*0 



2Asinh7t 



exp(7^) 



7 



2Asinh7* 



5 _i_ A 
1 + 



(4.9) 
(4.10) 



(4.11) 



(4.12) 



In Eqs. I|4.7II and l|4.11|) I u is the Bessel function of imaginary argument, I v (z) — {—i) v J v (iz) 3"i 
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V. DISCUSSION AND CONCLUSION 

Focussing on the expression (|4.12H for the first-passage-time distribution W(t) we note that the damping constant 
7 defines two distinct time regimes, I/7 setting the characteristic cross-over time. In the long time limit for jt 3> 1 
the damping constant controls the behavior of W(t). From Eq. i|4.12H we infer 

9 Axn +A/A /<y \ 3/2+A/2A 

W ® * F(1/2:A/2A) (1) -P^ 1 + A /A)*], (5-1) 

i.e., falls off exponentially with an effective damping constant 7(1 + A/ A) renormalized by the ratio A/ A of the 

nonlinear strength to the noise strength. We note that for A — ► the result is in accordance with the weak noise 
phase space derivation in Sec. IIIII In the intermediate time regime for 'yt <C 1 the damping constant 7 drops out and 
we obtain 

_ 2 OA/ 



W{t) * 1WA72A) -P(-o/2A*) ^ J ■ (5-2) 

For 2 At 3> the distribution VK(t) exhibits a power law behavior with the same exponent 3/2 + A/2 A as in the 
undamped case for 7 = 0. For weak noise this result is again in agreement with the estimate in Sec. IIIII 

In the short time limit W(t) vanishes exponentially and shows a maximum about the finite-time-singularity. In 
Fig. |3 we have depicted the first-passage-time-distribution as a function oft. In Fig. Q] we illustrate the behavior of 
W(t) in a log- log representation. 

In this paper we have extended the model discussed in ref. jl.'il to include a linear damping term. Not surprisingly, 
the damping changes the long time behavior of the physically relevant first-passage-time distribution. The finite- 
time-singularity occurring at time to in the noiseless case is still effectively resolved by the noise, becoming a random 
event, but the power law scaling behavior with scaling exponent a = 3/2 + A/2 A is limited to early times compared 
with the cross-over time I/7 set by the damping constant. In the long time limit beyond I/7 the damping gives rise 
to an exponential fall-off and the scaling property ceases to be valid. To the extent that the present simple model 
might apply to physical phenomena where damping is always present, we must conclude that an eventual power law 
scaling presumably is confined to a time window determined by the size of the damping. 
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APPENDIX A: EXACT SOLUTION OF THE FOKKER-PLANCK EQUATION 

In this appendix we discuss the exact solution of the Fokker-Planck equation in more detail. Denoting the nor- 
malized eigenfunction of H in Eq. (|4.4|l and the associated eigenvalues by \& n and A 2 E n /2, respectively, we obtain, 
incorporating the initial condition P(x,0) = S(x — xq) and the gauge transformation, the following expression for the 
distribution 

P(x, f)=£ e- AE ^ 2 e-^ 2 -^ 2A (x/x )- x ^ n (x)n(xo)- (Al) 

n 

By means of the transformation ^>(x) = X 1+X / 2A exp — r )x 2 /2AG(^/x 2 /A) it follows that G is a solution of the degen- 
erate hypergeometric equation [35L l36| . For the discrete spectrum we choose the polynomial form and further analysis 
shows that the eigenfunctions \& n are given in terms of the Laguerre polynomials [3^, |3(j . For the normalized 
eigenfunctions we thus obtain 



7 s3/2+A/2A r(n+l) 



1/2 

;c l + A/2A e - 7 ,V2A i l/2+A/2A (7a ,2 /A)j (A2) 



-AJ T(3/2 + A/2A + n)_ 

with discrete eigenvalue spectrum 

E n = 4n 7 /A + 2( 7 /A)(l + A/A) . (A3) 
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Inserting Eqs. (|A2|) and ((A3(l in Eq. (| Al|l and using the identity |2a, |3 



E 



T{n + a + 1) 



z n LZ(x)LZ(y) 



(sy«0 -/» e -, (8+y ) (1 -,) Ja(2 (^)i/2 /(1 _ z)) ; 



(A4) 



we finally obtain Eq. (|4. 1 1|) for P(x,t) and by the same analysis as in Ref. 13] the expression 1)4.12(1 for W(t). 
We note that for A — * 0, using Ii/ 2 {x) — {2/ttx) 1 ^ 2 sinh a: [33, the expression (|4.11|l takes the form 



P(x,t) 



_ttA(1 



-2-yt) 



- 1/2 






exp ^ 



7 (x — x e 7 *) 2 
A 1-e- 2 ^* 



exp 



7 (x + xoe- 7 *) 2 
A l-e" 2 T* 



(A5) 



i.e., the mirror case of the noise driven overdamped oscillator. Note 1) that for A — > the variable x lies on the 
noiseless orbit x — > xq exp(— 7 i) and 2) P vanishes for a; = 0. 



APPENDIX B: SMALL NOISE LIMIT - SADDLE POINT ANALYSIS 



In this appendix we perform for completion a weak noise saddle point analysis of the exact expression in Eq. H4.11|) 
along the same lines as in ref. .13']. This analysis requires that we consider both large order and argument of the 
Bessel function I v {x). This is easily done by Laplace's method using a convenient spectral representation mill 



i> + i/2)r(i/2) 

Inserting I(B1() in (|4.11|) we have 



f cosh(x cos 6) sin 2l/ 6 dO. (Bl) 
J o 



p( x t ) = 1 1 A C A/2A f l x \ e -(X/4A)(£ J +£i)/(xx sinhu) p jt/2 

4ttV2x A \\smh~ft; 



~ 2 \ 1/2+A/2A 



dQ S ' n _ Z' e (A/A)(logsin8+(l/2) cos 8/ sinh u) _|_ g (A/A)(logsin 8-(l/2) cos 8/ sinh «)\ (B2) 

sinh it V / 



where u is defined by sinhu = Asinli7t/(27xio), x = xexp(jt/2), and xq = xq exp(—jt/2). Setting f±(6) — logsin#± 
(l/2)cos0/ sinhw the saddle points for small A are given by cos(9± = ±exp(— u) for x > and cos8± = =Fexp(u) for 
x < 0. For x > we have f+(0+) = /-(<?-) = (l/2)(log(l - e~ 2u ) + e u / sinhu) and /£(0+) = /"(6L) = -cothu, and 
we obtain the weak noise result for x > 

/ ^ \ 1 / 2 e 7t/2 /-yf2^ _ e _2u )\ 1 / 2+A / 2A e -(A/2A)(i 2 +£^~2x£ cosh «)/iS sinh « 

^ ' ^ \47rA/ \ Asinli7i / (sinh u cosh u) 1 / 2 ^ 



For A — > the factor i +ig — 2iio coshu in the exponent in l(B3(l locks onto zero, thus setting x +Xq — 2xxq coshu = 
and inserting sinhu = Asinh( 7 i)/2 7 :r;Eo we obtain 

x 2 +xl - {{2xx ) 2 + (Asinh7t/7) 2 ) 1/2 = 0. (B4) 

Finally, setting x, = irexp(7i/2) and xq = xq exjp(—jt/2) we obtain after some reduction 

x(t) = ^Jx 2 cM-^t) - (W - cxp(-2 7 i)), (B5) 
which by simple inspection is equivalent to to Eq. 1(2.31 For 7^0 we have x = \J x^ — At; for A — > 0, x = xoe^' yt . 
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FIG. 1: In a) we show the time evolution of the single degree of freedom x. At times shorter than the cross-over time I/7 the 
variable x falls off exponentially. At times beyond I/7 the variable x reaches the absorbing state x — at a finite time to. In 
b) we depict the free energy F(x) driving the equation. For A = the free energy forms a confining harmonic well, for 7 = 
we have the absorbing state case discussed in ref. [l3|. In the general case the absorbing state x = corresponds to the sink 
in F(x). 




x x* x E=0 

FIG. 2: We show the topology of phase space. The bold lines indicate the zero energy sub-manifolds. The invariant heteroclinic 
orbits p = 273; and p = X/x passing through the saddle point FP at (x* ,p*) = [(A/27) 1/2, (27A) 1 / 2 ] have energy — 7A. The 
arrows indicate the direction of the flow. The long time orbit from xo to x passes close to the fixed point. The part of the orbit 
following the invariant manifold p — X/x and entering in our long time estimate is denoted I; the part of the orbit close to the 
p = 27a; manifold is denoted II. 
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fW(t) 




1/y 

FIG. 3: We sketch the first-passage-time distribution W(t) as a function of t. In the limit — > tO W(t) vanishes exponentially; 
about the finite-time-singularity W(t) exhibits a maximum. At intermediate times for "ft <JC 1 the distribution exhibits a power 
law behavior with scaling exponent 3/2 + A/2A. In the long time limit for ■yt S> 1 an exponential fall-off with time constant 
7(1 + A/A) characterizes the behavior of W(t) 




FIG. 4: In this figure we sketch the behavior of W(t) in a log- log plot. At intermediate times earlier that I/7 we have scaling 
behavior with exponent 3/2 + A/2A, corresponding to a constant negative slope. In the long time limit the curve dips down 
indicating the cross-over to exponential behavior. 



